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A large-n approximation to the S — 1 antiferromagnetic chain, using the symmetric tensor 
representation and its conjugate, is developed to order 1/n in order to calculate the magnon wave- 
function and to study the effect of modifying the exchange coupling from J to J' on a single link. 
It is shown that a magnon boundstate exists below the Haldane gap for arbitrarily small negative 
J' — J but only above a certain critical value of J' — J for positive values. In the former case the 
binding energy vanishes as (J — J') 
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I. INTRODUCTION AND REVIEW 



Many features of S = 1 Heisenberg antiferromagnetic chains have been understood in considerable detail using a 
combination of numerical techniques, approximate field theory methods and a model with bilinear plus biquadratic 
exchange for which the exact groundstate is a simple valence bond solid (VBS). jp There is considerable interest in 
the first excited state, the triplet magnon, separated from the groundstate by the Haldane gap. However, this state 
is not known exactly even for the VBS model. 

Investigations have been made of the Heisenberg model with the exchange coupling, J, modified to a different 
value, J', on a single link. It can be seen that a magnon boundstate exists, localized near the modified link, whose 
energy goes to zero in the limit J' /J << 1 or ,]'/,] » 1. However controversy remains over whether this boundstate 
continues to exist for | J' — J\ << J. The numerical results in Q suggest that a critical value of J' — J is required, for 
either sign of J' — J, in order for the boundstate to exist. On the other hand, the more detailed numerical analysis of 
H suggest that this is only true for J' > J while for J' < J a magnon boundstate exists for arbitrarily small values of 
J — J' . There is a well known theorem of one dimensional quantum mechanics which states that an arbitrarily weak 
attractive potential produces a boundstate. (For a general and rigorous version of this theorem see B.) Since the 
magnon behaves in many ways like an ordinary stable massive particle, one might expect that this theorem should 
be applicable to the modified link problem, but this connection has not been established. 
' The valence bond solid state was in fact encountered earlier |{| as the groundstate of the bilinear exchange model 
in a certain large-n limit. This is defined by placing the symmetric 2-tensor representation of SU(n) on the even 
sites and its conjugate representation on the odd sites. This is sometimes referred to as a type of "fermionic" large-n 
limit because it arises in a model with 2 "colors" of fermions and n "flavors" with a projection onto color singlets on 
each site with 2 particles on the even sites and 2 holes (n — 2 particles) on the odd sites. This must be distinguished 
from various other large-n limits including the fermionic case with n/2 fermion of each color on each site and the 
bosonic model with n bosons on even sites and n anti-bosons on odd sites. In || the model was solved only to leading 
order in 1/n. In this order, there remains a 3-fold groundstate degeneracy which is lifted in higher orders. The VBS 
groundstate is degenerate, to leading order in 1/n, with the two fully dimerized groundstates Js|. The Haldane gap, 
of 0(1) was calculated in j^j but a determination of the corresponding magnon wave-function requires going to order 
1/n, due to an infinite degeneracy at lowest order. 

The purpose of this paper is to develop this large-n approximation to next order in 1/n. We show that indeed the 
VBS state is the actual large-n groundstate, the dimerized states having energy higher by 0(l/n). [It was erroneously 
claimed in || that this splitting only occurs at (3(l/n 2 ).] We also solve for the magnon wave-function in the large-n 
approximation. In this approximation it corresponds to a type of soliton-antisoliton boundstate held together by a 
linear potential. Our results allow us to estimate the accuracy of the large-n approximation for the physical case 
n = 2. The results are somewhat encouraging. In any event, this approximation is useful because it qualitatively 
captures the essential physics. As an application, we show that the large-n approximation allows for a simple analytic 
treatment of the single modified link problem. We show that a boundstate exists for arbitrarily small negative J' — J 
but for positive J' — J only when J' — J is greater than a critical value, in agreement with Mallwitz and Wang 
with the binding energy vanishing as (J — J') 2 in the J' < J case. This result follows essentially from the theorem 
of one-dimensional quantum mechanics that an arbitrarily weak attractive potential always produces a boundstate. 
The simplicity and generality of the argument suggests that it may remain true for all n. 

For general n, the states on even sites are labeled by a pair of symmetric SU(n) particle indices: 
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|« >= |i< > 



(1.1) 



and on odd sites by a pair of symmetric hole indices: 

\ij >= \ji > ■ (1-2) 

For all n > 3 even and odd sites are inequivalent. However, in the n=2 case they are equivalent with the 3 states on 
each link corresponding to the usual S z eigenstates: 

I 11 > = |aa >= \S Z = 1 > 
T> =|i2 >HS z = 0> 

r> =| 11 >=|5 z = -i>. (i.3) 

For the 2-site model the Hamiltonian acts as: 

h\v m >= -(i/ n )[6ir, ml > +str, km > +5{r, ml > +sir, km >\. (i.4) 

For the case n = 2, 

H = Si - §2 — 1 (1.5) 

Throughout this paper we use the Einstein summation convention for a pair of repeated indices one upper and one 
lower. All SU(n) singlet states for the chain can be written with all upper indices contracted with lower indices. 
States with uncontracted upper and lower indices can be projected into irreducible representations of SU(n) by the 
usual processes of symmetrization and subtracting of traces. 

For the 2-site model it is quite easy to see that the groundstate is the SU(n) singlet > with energy -2 in the 



large-n limit. The factor of 1/n in (1.4) is canceled by 5\ = n. Now consider a chain of arbitrary length. We may 
represent an arbitrary SU (n) singlet state by drawing a line ( "valence bond" ) between pairs of sites (one even and one 
odd) representing each contracted pair of indices. In the large-n limit the energy of such a state is simply equal to (-1) 
times the number of nearest neighbor valence bonds. Hence there are three degenerate groundstates in the large-n 
limit, each of which has one nearest neighbor valence bond per link, on average. Two of these states are the dimerized 
states with a pair of valence bonds between each even site and the odd site to its right (or left). [See Figure [j].] The 
third degenerate groundstate is the valence bond solid state with a single valence bond between each neighboring pair 
of sites. [See Figure |J]. Only by going to higher order in 1/n will the degeneracy between the dimerized and VBS 
groundstates be lifted. We show in the next section that the VBS state is the true unique groundstate once 0{l/n) 
corrections are included. 

In the large- n limit, the lowest excited state has one fewer nearest neighbor bond and hence an excitation energy of 
1. Assuming the VBS wave- function at positive and negative spatial infinity, we see that all such states correspond to 
configurations with a string of dimers between the VBS states and dangling bonds separating the VBS regions from 
the dimer regions. [See Figure [|.] These dangling bonds correspond to solitons between the nearly degenerate VBS 
and dimer approximate groundstates. Such states are labeled by 2 indices, one upper and one lower corresponding to 
the two dangling bonds, >. Subtracting the trace, > — (Sj/n)\ k ,k > gives a state in the adjoint representation 
of SU(n). In the n = 2 case this reduces to the spin-1 (triplet) representation. To lowest order in 1/n, all such soliton- 
antisoliton (ss) states have the same excitation energy, +1. It is necessary to go to 0{l/n) to split the degeneracy. 
Note that, if the soliton and antisoliton are separated by a distance, x, then the energy will grow as x/n because there 
is a region of size x which is in the "wrong" (dimerized) groundstate, whose energy is higher by an amount of 0{l/n). 
Since the kinetic energy for the soliton and antisoliton is also 0(1 /n), it follows that the binding energy is 0(l/n). 
We essentially must solve a lattice version of the Schroedinger equation for a particle in a linear potential. In general 
this will give several soliton-antisoliton boundstates, each corresponding to an adjoint representation magnon. The 
number of stable boundstates grows with increasing n. 

In the next section we study the translationally invariant case, calculating the groundstate and magnon states to 
0(l/n). In Section III we consider the single modified link. A preliminary version of Section II appeared in B. 



II. GROUNDSTATE AND MAGNON STATES 
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A. Groundstate 



We label the candidate groundstates as follows. 



In -s. — I ''i'' 2 - - . . 

l u |-" J42«3 1 1*4*5 

I pn ^ _ | n«2 . . 13U . . 

\ ± ^^ <^ |'" U122 5 5*324 

IHP -V — I ^1*2 . . *3*4 . . 

I 1 ^-^ — |— 7*3*4 ! J*5*6 



> 
> 
> 



(2.1) 



Thus the \EO > and \OE > states are related by a translation of one site. \EO > is the state in which the double 
bonds have even sites on the left and odd sites on the right and vice versa for \OE > (see Figure Q). We expect the 
energy of the non-symmetric states to be identical and that of the symmetric state to be different (since no symmetry 
connects it to the non-symmetric states). Note that all three states have one bond per link on average. We now show 
that the VBS state |0 > is the true groundstate for large n 

We wish to calculate the energy of the symmetric groundstate including the first order correction. This is 



=< 0\H\0 > I < 0|0 > 



(2.2) 



Now, we first consider the action of the Hamiltonian on just two sites of the symmetric ground state. That is, we 
calculate H r \0 > where H r is defined to be the term in H acting on sites r and r + 1. We consider only the two sites 
on which H r acts. 



Hr 
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H r \ lll2 ;^ >= -(1 + 2/n)\ lll2 r* > -(l/n)dg| T 



> 



(2.3) 



(2.4) 



The first term is simply a multiple of the original state. The second term consists of a double bond and a delta 
function. This delta function serves to contract the two sites on either side of the two sites considered here thus 
producing a |/3 r _i > state (see Figure ^). Thus the action of the Hamiltonian on the entire symmetric ground state 
(L sites) is 



H\0 >= -L(l + 2/n)\0 > -{1/n) \Pr > 



(2.5) 



If we assume that the states |0 > and |/3 r > are orthogonal, 

< 0\H\0 >= -L(l + 2/n) < 0|0 > 

And so 

< 0|ff|0 > / < 0|0 >= -L(l + 2/n) 



(2.6) 



(2.7) 
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FIG. 3. a state 



In fact, the states |0 > and \/3 r > are not orthogonal. However, their overlap (divided by < 1 >) is order 1/n (see 
Appendix |A|) and thus the effect of the non-orthogonality is order 1/n 2 . 

We now consider the non-symmetric groundstate. As above, we wish to calculate 

E^o =< EO\H\EO > I < EO\EO > (2.8) 

which is the energy of the non-symmetric ground state including the first order correction. Now, in this case we need 
to consider the action of the Hamiltonian on two inequivalent types of pairs of sites, those possessing a double bond 
between each other and those not. For a chain of L sites, there are L/2 of each type (assuming L is even). 

First consider H r with r even. Then H r acts on a pair of sites possessing a double bond between each other. We 
need only consider the two sites on which the Hamiltonian acts. 

H r \ ili2 r* > = -(i/n)[sii\ mi2 r 2 > +si i 2 \ iim , mh > +< 2 u, iim > +^| lim ,* im >] 

= -(l/n)[n\ mi2 r 12 > +\hm, mil > +U i ; im >+n\ lim ; im >] (2.9) 

We see that all of the states in the above expression are simply the original state with a change of dummy variables. 
Therefore, 

H r \i lt2 ; 112 >= -2(1 + l/n)\ lll2 r* > (2.10) 

Now consider H r with r odd. Then H r acts on a pair of sites not possessing a double bond between each other. 
Again, we only consider the two sites on which the Hamiltonian acts. Remembering that the indices of the leftmost 
site are contracted to the left and those of the rightmost site are contracted to the right, 

HrWiJ 3U >= -(V™)[<^U, mi4 > +%Wm +<^U, i3m > +^|W 3m >] (2-11) 

In each of the states in the above expression, there exists a single bond between the two sites considered (that is, 
sites r and r + 1), a single bond between sites r — 1 and r resulting from one of the original dummy indices on site r 
remaining and a single bond between sites r + 1 and r + 2 resulting from one of the original dummy indices on site 
r + 1 remaining. Also, in each case, the delta function serves to contract the indices on sites r — 1 and r + 2. Thus, 
the action of the Hamiltonian on the state \EO > is 

H\EO >= -L(l + l/n)\EO > -(4/n) ^ \a 2r > (2.12) 

r 

where \a r > is defined in Figure [| If we assume the states \EO > and \ct2r > are orthogonal, 

< EO\H\EO >= -L(l + 1/n) < EO\EO > (2.13) 

And so 

< EO\H\EO > I < EO\EO >= -L(l + 1/n) (2.14) 

In fact, the states \EO > and \a r > are not orthogonal. However, their overlap (divided by < EO\EO >is order 1/n 
(see Appendix |X|) and thus the effect of the non-orthogonality is order 1/n 2 . Therefore, the above result holds to 
order 1/n. 



Comparing the values found in equations (2.7) and (2.14) we see that the symmetric ground state is the true ground 
state (having the lower energy) which agrees with known results. The energy difference per site is 1/n + 0(l/n 2 ). 
Considering states of the form 

> =|0>+(a/L)X>r > 

r 

|^2 > = \EO > +[b/(2L)] l«2r > (2.15) 

T" 

and calculating the ground state energies variationally produces the same results as those found above to 0(l/n) (see 
Appendix [b]) . Extrapolating the large-n approximation to n=2, and taking into account the shift of the Hamiltonian 
by a constant noted in ( |1.5| ), the estimates of the ground state per site to th and 1 st order in 1/n are and -1 
respectively. The numerical result for n — 2 is -1.401485. jq| Thus we see that the inclusion of the terms of first order 
in 1/n brings us closer to the numerical result. 
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FIG. 4. Soliton-antisoliton state. Note that an uncontracted upper index is represented by an upwards facing free link while 
an uncontracted lower index is represented by a downwards facing free link. 



B. The First Excited States 



Since our Hamiltonian counts the number of nearest neighbor bonds to leading order in 1/n, the excited states 
consist of an uncontracted lower index at some site and an uncontracted upper index at some other site (that is, 
they have one fewer nearest neighbor bond). These uncontracted indices can be though of as solitons/antisolitons 
interpolating between two "ground states" . We label such sites as in Figure [|. Note that r labels the uncontracted 
upper index and s labels the uncontracted lower index. We will assume that the indices on sites r and s are different 
to eliminate mixing with the singlet states (in which the indices on sites r and s are contracted). We define 

x = r — s (2-16) 

Note that r must be even and s must be odd implying that x is odd. Since the energy of the non-symmetric state 
is higher than that of the symmetric state by 1/n (+0(l/n 2 )), we expect the soliton-antisoliton pair to experience a 
linear (in x) confining potential like that encountered in quark confinement. 

Let us first consider H\f^s > for \s — r\ > 1. We need only consider the effects of H r —\, H r , -ff s _i, H s . The effects 
of the terms in the Hamiltonian acting on the other links have already been calculated. 

We will first consider the effect of H r -\. We let i\ be the index contracted between the sites r — 2 and i — 1, i% the 
index contracted between sites r — 1 and r, and ?3 the free index on site r. Now, 

H r -i\ lll2 r 3 > = -(i/n)[<^u™ 3 > +sz\nm m3 > +<tu, iam > +5? 2 \ nm ; 2m >] 
= -(i/n)[| mil r a > +n\ lim r s > +siiu 21 t2m >+k m : 3m >} 

= -(l/n)[(n + 2)\ mn r 3 > +<t U, iam >] (2.17) 

We recognize the first term as a constant times the original state (with renamed m) . The second term represents a 
double bond between states r — 1 and r while the delta function serves to create a free index (equal to the free index 
in the original state) on site r — 2. Thus, from this second term, we obtain the state \r — 2, s >. And so, 

H,._i|7Vs >= -(l + 2/n)|fTs > -(l/n)|r^2,s > (2.18) 

We now consider the effect of H r . Using the same indices as above and letting and i$ be the indices representing 
the double bond between sites r + 1 and r + 2, we calculate 

Hrl 1213 ^ >= -(l/n)[^n, mi5 > +5^ m , mt5 > +5%r i ^ iim >+5%\ i ^, iim >] (2.19) 

Now, in the second and fourth terms above, the index ii will contract to the left creating a bond, the index m 
represents a bond between sites r and r + 1, the indices i$ and (respectively) produce a bond between sites r + 1 
and r + 2, and the delta function serves to move the soliton (same free index) to site r + 2. Thus from these terms, we 
obtain the state \r + 2, s >. In the first and third terms above, the soliton (represented by the free index i^) remains 
on site r, the index m represents a bond between sites r and r + 1, and the indices is and ii (respectively) produce 
a bond between sites r — 1 and r + 2. Thus, in effect, we have replaced the bond between sites r — 1 and r and one 
of the bonds between sites r + 1 and r + 2 with one bond between sites r and r + 1 and one between sites r — 1 and 
r + 2. Thus, we have reduced the number of nearest neighbor bonds by one. This increases the energy by order 1. 
Therefore, we ignore these terms in our analysis. So, 

H r \f~s>=-(2/n)\r + 2,s> (2.20) 

By symmetry, 

H s \fTs > = -(1 + 2/n)\f%s > -(l/n)|r,s+ 2 > 
Hs-tlf^s > - -{2/n)\r,7^ 2 > (2.21) 
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Therefore, neglecting terms with energy higher by an amount of order 1, 

B\S?r > = -[(£ - 1)(1 + 2/n) - (|x| - l)/n]|f^ > ~(2/n)(|r + 2, s > + |r,7^ 2 >) 

-(l/n)(|r^2,s >+|r,s + 2>) (for |r - s| > 1) (2.22) 

We now consider the action of the Hamiltonian on the state \r, r + 1 > (that is, for x = 1). The action on |r, r — 1 > 
(x = —1) will follow by symmetry. We only need to calculate the effect of H r ; the rest are known. Assume the index 
i\ is contracted to the left, i± is contracted to the right, and %2, i% are free and consider 

HrU li2 ; 3U >= -(i/n)[sii\ mi2 r u > +5ii\ iim , mii > +<^u, i3m > +si% m , i3m >] (2.23) 

Now, the first term results in the state |r, r — 1 >. The second term results in a delta function times the symmetric 
ground state. The third term results in a state with a soliton on site r, an antisoliton on site r+1 and bonds between 
sites r and r+1 as well as r-1 and r+2 (this state has one fewer nearest neighbor bonds). The fourth term results in 
the state \r + 2, r + 1 >. Thus, ignoring states with higher energy of order 1 and assuming that ii ^ i^, 

H r \r,r + 1 >= -(l/n){\r,r~^ 1 > +\r + %r + l >) (2.24) 

and similarly 

H r \r,7^ 1 >= -(l/ri){\r - Xr - 1 > + |r,7+ 1 >) (2.25) 

At first glance, it seems as though our Hamiltonian is non-hermitian. However, we can see that this is simply due to 
the normalization of the \rys > states. The normalization for sites r + 1 through s— 1 is the same as that for an \EO > 
state of length |x| — 1. We see from Appendix [A] that this normalization is [2n(n + 1)] Cl^l i)/ 2 _ The normalization 
for the rest of the sites is the same as that for an |0 > state of length L — \x\ + 1 (where L is the total number of 
sites in the state |r, s >). It is shown in Appendix [A| that this normalization is ((n + l) i_ l T l+ 1 — l)/n. Thus the total 
normalization is 

< r7s|r7s >= [2n(n+l)]d a: l- 1 )/ 2 [((n + l) L -l x l +1 -l)/n] (2.26) 

In the large L and large n limits (taking L large first as usual), this reduces to 2" x > 1 '' 2 n L ~ 1 . Thus we define the 
"properly" normalized states 

\r, s >= 2-d a; l- 1 )/V 1 - L }/ 2 |7vs > (2.27) 

With this normalization, the Hamiltonian becomes 

{H + {L — 1)(1 + 2/n))\r, s > = [(|x| - l)/n]\r, s > -(V2/n)(|r + 2, s > +\r, s-2> 

+|r - 2, s > +\r, s + 2 >) for |x| > 1 
(fl"+(L-l)(l + 2/n))|r,r+l > = -(l/n)(\r, r - 1 > +|r + 2,r + l >) 

{H+{L-l){l + 2/n))\r,r-l> = -(l/n)(|r - 2, r - 1 > +\r, r + 1 >) (2.28) 

Now, we expect the first excited states to be translationally invariant. Therefore we define the translationally 
invariant states 



|x >= (1/V2L) V |2r,2r + x > (2.29) 



where the multiplicative constant ensures that < x|x >= 1. (This follows from the normalization of |r, s > since it is 
shown in Appendix |c| that < r', s'\rys > is of order n L ~ 2 or smaller if r' 7^ r and/or s' 7^ s.) 
Therefore, for |x| > 1, 

(H+(L~ 1)(1 + 2/n))|x > =(H + (L-1)(1 + 2/n))[(l/V5I) l 2r ' 2r + x >] 

r 

= (1/V2L) ^][((|x| - l)/n)|2r, 2r + x > -(V2/ri)(\2r + 2, 2r + x > 

r 

+\2r, 2r + x - 2 > +\2r - 2, 2r + x > +\2r, 2r + x + 2 >)] 
= [(|x| - l)/n]|x > -(2V2/n)(\x - 2 > +\x + 2 >) (2.30) 
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And for x = ±1 



{H+{L- 1)(1 + 2/n))|l > = -(2/n)| - 1 > -(2v^/n)|3 > 
(H+{L- 1)(1 + 2/n))| - 1 > = -(2/n)|l > -(2-\/2/n)| - 3 > (2.31) 

Although the assumption will prove to be a poor one in future calculations, if we assume that our wave function 
varies slowly, we can change our Hamiltonian into continuous form by Taylor expanding the states involved to obtain 

\x ± 2 >« \x > ±2(d/dx)\x > +2{d 2 /dx 2 )\x > (2.32) 

Therefore, 

\x + 2 > +|x - 2 >ra 2jac > +A(d 2 /dx 2 )\x > +... (2.33) 

Making this substitution into the expression for our Hamiltonian and ignoring the fact that the action of the Hamil- 
tonian is slightly different for x = ±1 than it is for \x\ > 1, 

H = -[(£- 1)(1 + 2/n) - 4\/2/n] + (|x| - l)/n- (8V2/n)d 2 /dx 2 (2.34) 

This is the Hamiltonian for a particle in a linear potential with x corresponding to the difference coordinate of the 
soliton-antisoliton pair. Thus, in this approximation, the soliton-antisoliton pair experiences a linear confinement 
potential similar to that encountered in quark confinement. If the pair becomes too separated, it splits into two 
pairs much like the quark-antiquark pair comprising a meson (this will be discussed in more detail later). The 
exact Hamiltonian (in the large-n limit) is of this form but the kinetic energy term is not equal to p 2 as the above 
approximation would suggest. 
We now define a general state 

|V; >= Yl ^( x )\ x > (2-35) 
x odd 

We then seek the lowest energy eigenstate of this form. This will be the first excited state. Thus, we wish to find the 
function ij)(x), defined on the odd integers such that 

(H-(E + 1 + 2/n))\ip > = SE\i> > 
=> i>( x )( H ~ ( E o + 1 + Vn))\x > =5EJ2 i>{x)\x > (2.36) 

xodd xodd 

where Eq — — L(l + 2/n) is the groundstate energy. Thus 

- (2/n)M-l)|l > +V(1)| - 1 >] - (2>/2/n)[^(l)|3 > +V(-1)| - 3 >] 

+ 4>{x)[{{\x\-l)/n)\x> -{2V2/n){\x + 2> +\x-2>)] = 5E^2^{x)\x> (2.37) 

xodd,\x\>l xodd 

We now take the overlap of both sides with the state < x'\. We have from the calculations of Appendix ^| that the 
overlap < x'\x > is order if x' — x and is of smaller order in n if x' ^ x. Thus, from the overlap with |1 > we 

obtain, to order l/n, 

- (2/n)</>(-l) - (2v / 2/n)V'(3) = 5Eip(l) V(3) = -(2ip(-l) + u5^(l))/(2v / 2) (2.38) 
And for |i| > 1 we obtain, to order l/n, 

- (2V2/n){ip{x - 2) + + 2)] + ((|a;| - l)/n)^{x) = 5E^(x) 

=> il>(x + 2) = [(\x\ - nSE - l)/(2V2)]ip(x) - xp(x - 2) (2.39) 

Now, since the Hamiltonian is an even function of x, we expect ip(x) to be either an even or an odd function of x. We 
can easily prove that the first excited state is even by a contradiction argument inspired by Feynman. Assume that 
the first excited state is odd (that is, tjj(x) is an odd function of x). Now define 
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9{x) 



—tp(x),x < 
ip(x), x > 



(2.40) 



which is obviously an even function of x. Now, the two states defined by these functions obviously have the same 
normalization. Thus, we can ignore this normalization when comparing the energies of the two states. The energy of 
the state \tp > relative to Eq + 1 + 2/n is 



= -(4/n)^(-l)^(l) + K(N _ l )l n )^ x f - C2V2/n)(^(x)^(x - 2) + ^(x)^(x + 2))] 

xodd, \x | > 1 

= (4/n)^(l) 2 + 2 ^ [((|a;| - l)/n)^(x) 2 - {2y/2/n){^{x)^{x - 2) + + 2))] 

xodd.x> 1 

while the energy of the state defined by g(x) is 

E g = -(4/n) ff (-l) S (l) + £ [((1*1 ^ l )l n )a( x f ~ (2V2/n)(g(x)g(x - 2) + g(x)g(x + 2))] 



(2.41) 



xodd, \x\ > 1 



-(4/n)V(l) ; 
- (8/n)V-(l) 2 



53 K(M " V/nMx) 2 - (2V2/n)(V(a:)V(x - 2) + VfrM* + 2))] 



a;od<i,a;>l 



(2.42) 



Thus the function defined by g(x) has a lower (or equal if = 0) energy than the state \tp >. This contradicts the 
claim that > is the first excited state. Thus the first excited state is an even function of x. 

Therefore, to find the first excite d sta te, we can set = 1 and need only determine ip(x) for positive x. Since ip 
is even, 1) = = 1 and so ( [2.38 ) becomes 



^(3) = -(ndE + 2)/2V2 



(2.43) 



Equation 2.43 gives 0(3) in terms of "0(1) = 1 an< i nSE. Equation [2.39 is a recursion relation defining the ip(x) for x 
odd. Now, in order for the state \tp > to be normalizable, we need the sum J2T=o ^{^ z + I) 2 to converge. We expect 
that this will only occur for discrete values of nSE. Thus, using the above recursion relation, we plot the value of a 
partial sum of the above series for various values of nSE. This plot is shown is Figure ||. 
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FIG. 5. Log of partial normalization sum versus nSE. The log of the sum is plotted to allow the wide range to be displayed 
more clearly. 
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FIG. 6. A blow up of Figure |E| near the first minimum. 



. 8 



H °- 6 



0.4 



0.2 



10 



15 



20 



X 

FIG. 7. Ground state wave function. Since the wave function is even, it is only plotted for positive x. Note that the wave 
function starts to diverge at about x = 19. This is due to the fact that the eigenvalue found is merely a numerical approximation 
to the true eigenvalue. 
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We find that the first excited state has nSE w —3.747 (see Figure [jj). The wave function for this energy is shown 
in Figure |^. Thus the energy gap of the first excited state to order 1/n is 

A = E 1 - E = 1 + 2/n + 5E w 1 + 2/n - 3.747/n = 1 - 1.747/n (2.44) 

Thus, for n=2, A w 0.1265. Numerically and experimentally, this value is found to be 0.41050 JjJ. Since we obtain 
A = 1 to leading order, and A w 0.1265 to order 1/n it seems plausible that were one to include higher orders, the 
calculated value of A might approach 0.41050 in an oscillatory manner. 

Using our numerical results to order 1/n we can also calculate how many bound states exist. We assume that a 
necessary and sufficient condition for the stability of a given state is that it has less than twice the energy of the first 
excited state (otherwise, it would decay into two or more states of lower energy). Thus, the m th excited state will be 
stable if and only if 

1 + 2/n + (nSE) m /n < 2(1 - 1.747/n) 

=> (n5E) m < n - 5.494 (2.45) 

Now, we find that the second excited state is a parity odd state with (nSE)2 ~ 0.467. Thus, the second excited state 
is stable if and only if 

0.467 < n- 5.494^ n> 5.96 (2.46) 
Thus, for n = 2 there exists only one stable magnon in agreement with cr-model, numerical and experimental results. 



III. SINGLE MODIFIED LINK 



We now consider the case in which one of the links in our chain has a coupling which may differ from that of the 
other links. In the following discussion, we will find it useful to use the coordinates 

x = s — r 

y=s+r-l (3.1) 

rather than r and s to label our magnon states. The fact that r must be even and s must be odd implies that either 
x — Ap + l,y = Aq or x = Ap + 3,y = Aq + 2 where p and q are both integers. Thus although x must be odd and y 
must be even, not all combinations of odd x and even y are allowed. 

We define the coupling between sites and 1 to be J' and between all the other sites to be J. As in the previous 
calculations, we will choose our units of energy so that J = 1. We assume that J' — 1 is 0(1 /n). This implies that 
the ground state will remain unchanged to leading order in 1 /n and that its energy will be (to order 1 /n) 

E = — (L + J' — 1)(1 + 2/n) (3.2) 

Note that this redefinition of Eq corresponds with our earlier definition in the case J' = 1. We now define Ho to be 
the Hamiltonian without the modified link with t he new generalized Eq substituted for the old (that is, we define 



H - E - (1 + 2/n) to be the right sides of (2J2 



We can write the true Hamiltonian (with the modified link) as Hq plus terms correcting for the modified link. That 
is, we can write the true Hamiltonian H as 

H = H + XH k + XV (3.3) 

where A = 1 — J' = 0(1/ n) and 

H k \r, s > = (V2/n)\r + 2, s > for r = 0, s/1 
= (V2/n)\r,s - 2 > for s = 1, r/0 
= (l/n)[\r + 2,s > +\r, s - 2 >] for r = 0, s = 1 
= otherwise 
V\r, s > = -(1 + 2/n)\r, s > , for r < 0, s > 1 
= +(1 + 2/n)\r, s > , for r > 2, s < -1 

= otherwise (3.4) 
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FIG. 8. Modified Potential V. An open circle denotes a coordinate at which V has the value +1 and a closed circle denotes 
a coordinate at which V has the value -1 



Thus XHk is a modified kinetic energy operator and AV is a modified potential. Now, since we are considering the 
large n limit, we can neglect Hk and drop the 2/n ter m in V since these are lower order in n. [Recall that A is 0(l/n).] 
Thus, if we now change coordinates to x and y, (3.4) becomes to leading order in 1/n 



H k \x,y > = 
V\x,y > = — \x, y > , for x > 1, |y| < x — 1 
= > i for x < -3, |y| < -(x 

= otherwise 



•3) 



(3.5) 



Thus, the modified link not only breaks translational symmetry, it also breaks parity symmetry (reflection about a 
site) since V is not symmetric with respect to parity. Our potential arises from the fact we have defined our energy 
such that zero energy corresponds to the modified link being outside the soliton-antisoliton pair where there is a single 
bond between neighboring sites. If the modified link is instead between the soliton and the antisoliton, it is either 
between two sites not bonded together (in this case, V=-l) or between two sites with a double bond (in this case 
V=+l). We see from (3.5) and Figure || that if the modified link is between the soliton and antisoliton (that is, V is 
non-zero) then, it must lie between two sites not bonded together if x > 1 and on a double bond if x < —3. 
Our goal is to determine the nature of possible bound states of the Hamiltonian given above. Let us define 



iVto >= X! ipo(x)\x,y> 
allowed x,y 



(3.6) 



11 



^ "jl. 



r=-4 1 

(a) 

S= . 3 ^ ^ ^ 

| 1 r=4 

(b) 

FIG. 9. (a) Typical configuration with r < 0, s > 1 (=>■ x > 1). There is no bond on the modified link 01. (b) Typical 
configuration with r > 2, s < —1 (=> x < —1). There is a double bond on the modified link 01. 



to be the translationally invariant state we found in the case of the unmodified Hamiltonian, normalized so that 

£ M*) 2 = 1 (3.7) 
x odd 

We see that in the limit of large L (recall L is the length of our chain), j^o > has energy Eq + A (note the change in 
the definition of Eq here) since the state is spread along the entire length of the chain (that is, it is independent of 
y) while the modified link occurs only between two sites and thus its effect becomes negligible. Therefore, we expect 
that an eigenstate of the modified Hamiltonian will be bound (that is, be normalizable) if and only if it has energy 
less than Eq + A. Let |0 > be such a bound eigenstate where 

|0 >= £ <Kx,y)\x,y> (3.8) 
allowed x,y 

Since |0 > is a bound state, we know that < 0|0 >< oo. Thus, the energy of the state |0 > is 

E^ = < <p\H \(j> > / < 0|0 > +A < 0|V|0 > I < 0|0 > 

> E + A + A < 0|V|0 > / < 0|0 > (3.9) 

where the inequality is due to the fact that, by definition, A is the lowest energy eigenstate of Hq — Eo (and therefore 
Eq + A is the lowest energy eigenstate of Hq). Now, we know that 0(cc,y) must decay rapidly at large due to the 
linear potential present in the Hamiltonian. We see from ( |3.5[ ) and Figure || that for small x, the potential V is only 
non-zero for small y. Since our perturbation is small (order A), we expect that for small x and y, ~ if) (that is, the 



corrections are order A and thus can be ignored to leading order). Therefore (3.9) becomes, to leading order in A, 



^>£ o + A + A<V>o|^|V>o>/<0|0> (3.10) 

The fact that ipa( x ) is independent of y makes < - 0o|^|'0o > particularly easy to calculate. Since tl>o(x) is an even 
function of x, both signs of x are weighted evenly in the sum < 7/)o|V r |V'o >• Also, the fact that ipo( x ) is independent 
of y means that the weighting is the same (for a given value of no matter where the modified link is. Now, for a 
given \x\, there is always one more link with no bond than with a double bond. This can be seen easily in Figure]^. 
Thus, for each |x|, the sum over +x and — x and the respective allowed values of y in < ■0o|T^|'0o > equals -1 since V 
takes on the value -1 one more time than it takes on the value +1. Thus, 

< 4> \V\tJj q >= - £ ^ (x) 2 =-l/2 (3.H) 

a:>l Odd 



since this sum is equal to half the normalization sum in ( p.7[ ). Therefore, 

E^ > E Q + A - A/(2 < 0|0 >) (3.12) 

Since < 0|0 >> 0, if A < then E$ > Eq + A which is a contradiction since this implies |0 > is not a bound state. 
Thus, in the small A limit, it is only possible for a bound state to occur if A > (that is, if J' < J). However, we do 
not yet know that a bound state does exist in this case. To see this, we will use a variational argument. 
We choose as our variational state the normalized wavefunction, 

\ip> = ^2 H x >v)\ x >y > 

allowed x,y 

where ip(x, y) = Aip (x)e~ eM and A -> Vie as e — ► (3.13) 
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To see how the normalization constant A (chosen so that < ip\i/j >= 1) was obtained, recall that either x = 4p+ 1, y — 
Aq or x = Ap + 3, y = Aq + 2 where p and q are integers. Thus, 



allowed y 



geZ 9=0 



oo 

Oi? =2e- 4e Xe- 8£|<?l =2e- 4 7(l-e- 8e )«l/(4e) (3.14) 

9=0 



which gives us our normalization constant A in the small e limit. We choose the trial wave function given above since 
it reduces to the translationally invariant wave function in the limit e — > (we expect e — ► as A — > 0). We will 
minimize the energy of this trial wave function with respect to e. 

In terms of the coordinates x and y and using the same normalization for the \x,y > states as was used for the 
|r, s > states, Ho is given by 

{H -(E + l + 2/n))\x, y > = ((\x\ - l)/n)\x, y > -(V2/n)(\x + 2, y + 2 > +\x + 2, y - 2 > 

+ \x-2,y + 2 > +\x-2,y-2>) for |x| > 1 

(Ho -(Eo + l + 2/n))|l, y > = -(V2/n)(\3, y + 2> +|3, y-2>) 

-(l/n)(|-l,y + 2>+|-l,y-2>) 
(H -(E + l + 2/n))| - 1, y > = -(V2/n)(\ - 3, y + 2 > +| - 3, y - 2 >) 

-(l/n)(|l,y + 2>+|l,y-2>) (3.15) 

We would like to separate H into its a; and y components. One possible method for doing this would be to introduce 
lattice derivatives. However, a problem arises when we atte mpt this. We would like to add and subtract the states 
\x + 2, y > and \x — 2, y > from the right hand side of (|3. 15 ) . However, given that the state \x,y > exists, neither of 
the states \x + 2, y > or \x — 2, y > can exist (this follows from the permissible values of x and y). However, the states 
I — (x + 2),y > and | — (x — 2), y > do exist. Although adding and subtracting these states from the right hand side 
of (3.15) may seem counterintuitive, the fact that tp(x, y) is an even function of x allows us to effectively separate Ho 
into its x and y components in this way. Thus, we rewrite the \x\ > 1 portion of (3.15) as 

(Ho - (Eq + 1 + 2/n))\x, y > = ((\x\ - l)/n)\x, y > -(V2/n)(\x + 2, y + 2 > +\x + 2, y - 2 > -2\ - (x + 2), y > 

+ \x - 2, y + 2 > +\x - 2, y - 2 >) - 2| - (x - 2),y >) 
-(2V2/n)(\ - (x + 2), y > +| - (x - 2), y > -2\x, y >) - (*V2/n)\x, y > 
= ((\x\ - l)/n - AV2/n)\x, y > -(V2/n)(A y \x + 2,y> +A y \x -2,y>) 

-(2V2/n)A x \x iy > (3.16) 

Where 

A x \x,y > = \-(x + 2),y>+\ - (x-2),y> -2\x,y> 

A y \x,y> = \x,y + 2 > +\x,y- 2 > -2| - x,y > (3.17) 

are our lattice second derivative operators. Similarly, 

(H -(Eo + l + 2/n))|l, y > = -(V2/n)A y \3, y > -(l/n)A y \ -l,y> -(2y/2/n)\ - 3,y > -(2/n)\l, y > 
(H Q -(E a + l + 2/n))\ - l,y > = -(V2/n)A y \ - 3,y > -(l/n)A y \l, y > -(2y/2/n)\3, y > -(2/n)\ -l,y> (3.18) 

We can now roughly separate Hq into its x and y components. Define Hq as the sum of the terms in ( |3.16 ) and 
( |3.18 ) involving A y and Hq as the sum of the remaining terms. Since ipo(x) is an even function of x and is independent 
of y, A y \ipo >= 0. Thus H^tpo >= and so Hq is precisely the Hamiltonian we encountered in the case J' = J up 
to a constant relating to the redefinition of Eo- 

We now calculate the expectation energy of our trial wavefunction. 
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E c = < i/>\H\ip > 

= E + A+ < ip\H%\ilj > +\ < ip\V\ip > (3.19) 

Now, although A y is defined as an operator on the state \x,y >, through a change of variables, we can treat it as an 
operator on the coefficients ip(x,y). For y > (note that y is always even), 

A y tp{x,y) = VteMx)[e~ <y+2) + e~< y -V - 2e~ ey ] 
= iJKx,y){e~ 2e +e 2e -2} 

= V(x,2/)[4e 2 + 0(e 4 )] (3.20) 
where we have used the fact that ipo{x) is an even function. Similarly, for y < 0, 

Ay1>{x,y) = ^(x 1 y)[4e 2 +0(e i )} (3.21) 

Now, for y = 0, 

A y ^(x,0) = vTeMx)[2e~ 2e - 2] 

= ^(x,0)[-4e + 4e 2 + O(e 3 )] (3.22) 

Notice that the case y = agrees with the case y ^ in the e 2 term. Therefore, dropping terms of order e 3 and higher 
and using the fact that i/j{x, y) is an even function of x we see that 

< > = ~(V2/n)4e 2 £ [^(x, y)i>(-(x - 2), y) + ^{x, y)^(-{x + 2), y)] 

allowed \x\>i,y 

-(V2/n)(-4e) £ [^(x,0)V(-(x - 2),0) + i/>(x,0)il>(-(x + 2),0)] 
allowed \x\>i 

-4e 2 J2 [(V2/n)V(l,y)V(-3,y) + (l/n)^(l,y)^(l,y)] 

allowed y 

+4e[(V2/n)^(l, 0)V(-3, 0) + (l/n)V(l, 0)^(1, 0)] 

-4e 2 Yl [(\ / 2/n)V(-l,y)V'(3,y) + (l/n)V(-l,2/)V(-l,2/)] (3-23) 
allowed y 

Now, we can use the fact that x, y) = ip(x,y) to convert sums over all x to (two times) sums over positive x. 
Also, the x — ±1 cases missing in the first two sums above, are present in the third and fourth sums. Finally, in the 
second sum of (|U|), we must remember that only certain values of x are allowed if y = 0. The allowed values are 
...,-7,-3,1,5,9,.... Thus, using the fact that i^o(x) is an even function of x, we can convert a sum over these allowed 
values into a sum over all odd positive x. Thus, we may simplify the above expression to: 

< m V oH > = "(V2/n)4e 2 [4 £ ^(x, y)^{x + 2, y)] 

allowed x>l,y 

+(V2/n)(4e)[2 ^ ^(s, 0)i/>(x + 2, 0)] 
odd x>i 

-(l/n)4e 2 [2 £ ^(1,^(1, »)] 
allowed y 

+(4e/n)V>(l,0)V(l,0) (3.24) 

Now, a factor of e occurs from the normalization of ip which is canceled when a sum over y is performed. Thus, all 
the terms in the above end up being 0(e 2 ). Thus, the above becomes 

< ii\H v a \ij > = -{V2/n)16e 2 ^ ip {x)ip {x + 2) + (V2/n)8e ^ AeipoWM* + 2 ) 

odd x>i odd x>i 

-(8e 2 /n)Ml) 2 + (16e 2 /n)^ (l) 2 
= (e 2 /n)[16V26 + 8^o(l) 2 ] (3.25) 
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where 



9 = Mx)M* + 2) > (3.26) 

odd x>i 

It is shown in Appendix [d| that ipa is everywhere positive and so we know 6 is positive. 

We now consider the term V, Using the same reasoning as when we were calculating < 0|V|0 > (that is, for small 
x, y and e, ijj(x,y) ~ Aipo(x)), we see that 

< ip\V\ip >« A 2 < ip \V\ip >= (4e)(-l/2) = -2e (3.27) 

Therefore, the expectation energy of our trial wave function is 

E € = E + A + 8Ke 2 /n - 2Ae 
Where k = 1\p2B + "0o(l) 2 > (3.28) 

In order to minimize this quantity with respect to e, we set dE/de = to obtain 

= 16«e/n- 2A 

=*> e = tiA/(8k) (3.29) 

which gives an energy of (to leading order in n) 

E mm = E Q + A - n\ 2 /{8 K ) <E a + A (3.30) 

Thus, we know that the true magnon state of the modified link Hamiltonian has an energy which is less than Eq + A 
and thus is bound. Note that our minimizing value of e is positive which is necessary for our trial wave function to be 
normalizable. Numerically, 9 is calculated to be approximately 0.279 and thus k is approximately 1.134. Therefore, 
for n = 2 and A > 0, 

E mm ^E + A~ 0.221A 2 (3.31) 



So far we have presented the wave-function of (3.13) as merely a variational one so that the energy E m i„ is just 



an upper bound. However, we expect the variational wave-function to become sufficiently accurate that the actual 



boundstate energy is given by (3.3C) in the limit A — > + . This follows from the following observations. First of all, 
ip(x, y) is an exact eigenfunction in the translationally invariant case, A = 0, for any e. Secondly, even for A ^ 0, ip( x i y) 
satifies exactly the Schroedinger equation outside the dotted lines in Figure ||. For any fixed value of x and sufficiently 



small A [and hence e as determined by ( p. 29 )] ijj(x,y) is nearly constant as a function of y between the dotted lines 



of Figure^. The true boundstate wave function must also have this property in order for < i(]\Hq\iI> > to be small. 
Thus the actual value of the wave-function in this region is not important as long as it is nearly constant and joins 
smoothly with the wave-function outside the dotted lines, properties enjoyed by rp(x,y). At \x\ >> 1, ip(x,y) goes to 



rapidly. [From the continuum form of the Hamiltonian in (2.34) we may estimate ipo{%) e x 1 ^ 3 ' 2 1 Again 



we expect the actual boundstate wave-function to have this property so ip(x, y) should be a sufficient approximation 



in the large \x\ region. Therefore we expect ip(x,y) [with e given by (3.29)1 to be a good first approximation to the 



actual boundstate wave- function at small A and the energy of (3.3C) to be asymptotically correct. The situation is 
similar to the case of one-dimensional quantum mechanics with a weak potential, the single dimension corresponding 
to y. Despite the fact that -0 depends on 2 co-ordinates x and y, the present problem is rather different than the 
two-dimensional case with a short-range potential due to the confining x-dependent potential. Essentially, x acts as 
an internal degree of freedom of the magnon whereas y represents its location. In particular, the quadratic dependence 
of the binding energy on the strength of the potential also occurs in the one-dimensional case. 

Thus, if we assume that E m i„ is not merely an upper bound on the energy but actually a good approximation to 
it, we can plot the energy as a function of A = (J' — J)/ 'J. This is done in Figure [l0| Note that the critical value of 
A > at which the energy becomes negative (implying the existence of a bound state) has not been calculated and so 
no significance should be attributed to the value indicated on the plot. Also, the exact form of the energy dependence 
for A > is unknown. In Figure [w] it is plotted with the same quadratic dependence as in the case A < 0. Based on 
our experience in the previous section with extrapolating large-n results to n=2, we might expect the pre-factor of 



.221 in Eq. (3.31) to only be accurate to within a factor of 2 or so. However, the fact that the functional dependence 



is quadratic in (J' — J) is likely to be exact, since it is a general result for one-dimensional quantum systems, as 
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(J'-J)/J 




FIG. 10. E m i„ versus (J' — J)/ J. The critical value of J' > J for which the energy becomes negative has been assigned an 
arbitrary value. The form of the energy dependence for J' larger than this critical value has been assumed to be the same as 
the dependence for J' < J. 

emphasized above. The numerical results in || are probably consistent with quadratic behaviour at small negative 
(J' — J) but with a smaller pre-factor. 

Therefore, we have shown that in the case of a modified link, in the limit J' — > J, a bound state exists if and only 
if J' < J. This results from the fact that between the antisoliton and soliton there is one more pair of neighboring 
sites with no bond between them than with a double bond. 
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APPENDIX A: CALCULATION OF OVERLAPS ARISING IN DETERMINATION OF TRUE GROUND 

STATE 

1. Calculation of < 0|0 > 

For the symmetric ground state, we first consider just two of the sites. 

= (n + 2)#a£+$$ (Al) 

We then consider three sites. 

= [(n + l)(n + 2) + l]^ + ^ (A2) 
We see that this pattern will continue. That is, for L sites, with L even, we obtain 

<u\o>=*l%€:i+s? l+ Ai +i ( A3 ) 

Where is given by the recursion relation 



1G 



a L+1 = (n + l)a L + 1 (A4) 

However, this is simply a geometric series. 

ai = 1 

a 2 = (n + 1) + 1 

a 3 = (n + l) 2 + (n+ 1) + 1 

L-l 

aL = 5Z(n+l) J = [(n + l) L -l]/n (A5) 



i=0 



Thus for a finite chain of L sites, 



< 0|0 >= [((n + 1) L - l)/n]^<5^ + ( A 6) 

If we now wish to impose periodic boundary conditions, this is equivalent to replacing by i\ and jl+i by ji in 
the above expression. Doing this, and contracting the indices yields 

< 0|0 >= [((n + 1) L - 1) /n]n + n 2 = (n + l) L + n 2 - 1 (A7) 

Taking L — > oo first, then large n, this may be approximated as: 

<0|0>=(n+l) L . (A8) 



2. Calculation of < EO\EO > 

To calculate this matrix element, we can use the previous result. Each double bond can be thought of as a two site 
|0 > state with periodic boundary conditions. Thus, for L even 

< EO\EO >= \{n + I) 2 + n 2 - 1] L ' 2 = \2n{n + 1)] L/2 (A9) 



3. Calculation of < a r \a r > 

We again use our previous results to calculate this matrix element. The four sites which comprise the soliton- 
antisoliton pair can be thought of as a four site |0 > state with periodic boundary conditions while the remainder of 
the \a r > state is identical to the \EO > state. Thus, for L even 

< a r \a r >= [(n + l) 4 + n 2 - l][2n(n + l)]^" 4 )/ 2 (A10) 



4. Calculation of < /3 r \f3 r > 

The overlap of \(3 r > with itself consists of the overlap of a two site \EO > state (or equivalently a two site |0 > 
state) with itself and an (L — 2) site |0 > state with itself. Thus, for L even 

< (3 r \(3 r >= 2n(n + l)[(n + l) i_2 + n 2 - 1] (All) 
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|a r oO= ... ^ ^ s , r\ ^ ... ^ ^ , . % ^ ^ ... 

I r V r s 

FIG. 11. Double alpha state 

5. Calculation of < Q\0 r > 

We need only consider the four sites which comprise the soliton-antisoliton pair. The remainder of the matrix 
element is identical to < 1 >. 



H2l3 I )»4»5 \3lj2 ' 



(ft&sl' + + sfsfsi 5 + <)■;■&!■■■ &';■) 

v 33 34 *5 J3 «5 .34 J3 «B 34 l 5 ' 

= (2n 2 + 6n + 4)<J£<$|? + 2(n + l)^^ 

= 2(n + l)[(n + 2)^+^^] (A12) 

We now note that this is precisely 2(n + 1) times the expression we obtain from a two site < 0|0 > aside from a 
renaming of indices. Thus, since the remaining L — 4 sites of the \(3 r > state are identical to those for the |0 > state 

< 0|/3 r >= 2(n + l)[(n + 1) L ~ 2 + n 2 - 1] « [2/(n + 1)] < 0|0 > (A13) 



6. Calculation of < EO\a r > 

As in the previous calculation, we need only consider the four sites which comprise the soliton-antisoliton pair. The 
remainder of the matrix element is identical to < EO\EO >. 

< il<2 ,ni 2 r i ,i 3 i 4 \nh, hh m34 > = + + ttji^A + + S X) 

= [(2n + 2)$][(2n + 2)*£] 

= 4n(n + l) 2 (A14) 

Recall that for a four site \EO > chain, 

< EO\EO >= An 2 {n + if (A15) 
Thus, adding in the other L — 4 sites, we obtain 

< EO\a r >= (l/n) < EO\EO >= (l/n)[2n(n + 1)] L/2 « (l/n) < EO\EO > (A16) 



7. Calculation of < a r |a s > and < EO\a r a s > for r — s\ > 3 

In this case, the two soliton-antisoliton pairs do not overlap. Thus, comparison with the calculation of < EO\a r > 
yields 

< a r \a s >= (l/n 2 )[2n(n + 1)} L / 2 « (l/n) 2 < EO\EO > (A17) 



8. Calculation of < at\a r a s > for |r — s| > 3, \r — 1\ > 3, Is — 1\ > 3 



Here |a r a s > is defined in the obvious way in Figure 11. In this case, none of the soliton-antisoliton pairs overlap. 
Thus, comparison with the calculation of < EO\a r > yields 



< a t \a r a s >= (l/n 3 )[2n(n + 1)} L/2 



(A18) 
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FIG. 12. Double beta state 



9. Calculation of < /3 r |/3 s > and < 0\/3 r (] s > for \r — s\ > 3 



As in section A 7, comparison with the calculation of < 0|/3 r > yields 

< &.|& >= 4(n + l) 2 [(n + l) L - i + n 2 - 1] w [2/(n + l)] 2 < 0|0 > (A19) 

10. Calculation of < /3 t |/3 r /3 s > for |r - s| > 3, \r - 1\ > 3, |s - t| > 3 

Here |/3 r /3 s > is defined in the obvious way in Figure [l2| Comparison with the calculation of < 0\(3 r > yields 

< Pt\PrPs >= 8(n + l) 3 [(n + 1) L ~ 6 + n 2 - 1] (A20) 

APPENDIX B: VARIATIONAL VERIFICATION OF TRUE GROUND STATE 

We consider states of the form 

\tpi > = |0 > +a\0 > 

|^2 > = | SO > > (Bl) 

Where 

|/3 > =(1/L)53|A.> 

r 

|a> =[l/(2L)]^|a ar > (B2) 

r 

We will see that variations in the parameters a and b do not affect the energies of the states at order 1 /n thus verifying 
that the symmetric ground state is the true ground state (in the large n limit). Throughout this section, we will 
make the assumption that for large L, when the Hamiltonian acts on \a r > or \/3 r >, we can neglect the terms of the 
Hamiltonian acting on the five sites r — 1 to r + 3 since their effect is suppressed by a factor of L compared to the 
other sites (of which there are L — 5). 

For the symmetric ground state, we wish to calculate 

£i=<^i|JT|^>/<^i|Vi> (B3) 
We already know H\0 >. Now, for |r — s| > 3, 

H s \f3 r >= -(1 + 2/n)|/3 r > -(l/n)|A-ft > (B4) 
Where \(3 r f3 s > is defined in Appendix Thus, using our large L assumption, 

H\0r >= -L(l + 2/n)\0r > -(1/n) ^ > (B5) 



Thus, we see that the normalization conditions ((A7), (A13) and (A19)) and the action of H are consistent with 



\f3 [2/(n+ 1)]|0 > (B6) 

Although this equation is not literally true, the variational calculations give the same result as if it were true. 
Therefore, | Vx > just acts as a scalar multiple of |0 > and so 
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* ] 

~ s \- 11 iu ^4 ^ ... ^ 

/ — t u r 



r 



j 

t,u 



1 ^ ... ^ " ' h ^Jl 



t | r s 

FIG. 13. Site and index labeling when t < u < r < s. 

< > / < <Ai|V>i >=< 0|#|0 > / < 0|0 >= -L(l + 2/n) (B7) 

to first order in 1/n. 

For the non-symmetric ground state, we need to calculate 

E 2 =< t[>2\H\iP 2 > I < Hfo > (B8) 

Calculations analogous to those for the symmetric yield groundstate demonstrate that the normalization conditions 
((|A|), ( |Aie| ) and ( |A17| )) and the action of H are consistent with 



| a >« (l/n)\EO > (B9) 

Again, this equation is not literally true; the variational calculations just give the same result as if it were true. 
Therefore, \ip2 > just acts as a scalar multiple of \EO > and so 

< ip2\H\i/j 2 > j < ip 2 \ip2 >=< EO\H\EO > j < EO\EO >= -L(l + 1/n) (BIO) 

to first order in 1/n. 

Thus, variations in a and b do not affect the energies of the states | -01 > an d |*02 > a t order 1/n as claimed. 

APPENDIX C: CALCULATION OF OVERLAPS FOR SOLITON-ANTISOLITON STATES 

In this section, we will only consider states in which the soliton is to the left of the antisoliton (that is, \i\s > for 
which s > r). The other cases are easily seen to produce analogous results. In the forthcoming calculations, we will 
find the overlap of 2-site wave functions, one corresponding to a single bond and the other to a double bond, to be 
useful. This is just 

< ni2 ,m 2 \nn, J2j3 >= + + = 2 ( n + 1 H ( C1 ) 

This result can obviously be extended to longer chains. 



< 



U«2 lx-2lx-l I 3233 

nil2 ■■■ !Jaj-2ii!-l \3l32 > ■•■3x-23a 



^3.-13. >= [ 2 („+ (02) 



Where we note that in the above, x must be odd. 

Let us first calculate the overlap < fys\t, u > when t < u < r < s. We label indices as in Figure |l^. We define 
x' = u — t, x — s — r and d = r — u + 1. From ( |C2| ), the overlap for sites t + 1 through it — 1 is 

[2(n + l)]( x '- x )/ 2 4" (C3) 

Similarly, the overlap for sites r + 1 through s — 1 is 

[2(n + l)]( x - x )/ 2 ^ (04) 

From Appendix ^ we know that the overlap for sites u through r, 

[(^ + l) d -l)/n}Si : SX+Sj : SZ (05) 

and the overlap for sites s through t (periodic boundary conditions) is 
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FIG. 14. Site and index labeling when t < r < u < s. 
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FIG. 15. Site and index labeling when t < r < s < u. 



Ju 



[((" + l) L - d - x - x ' +2 l)/n]5ll5l + 8pl (C6) 

Therefore, the entire overlap is 

[2(n + l)](^'- 2 )/ 2 ^<5j : [(((n + l) d - l)/n)<t<^ + #£][(((„ + i)*-<*-*-» '+ 2 - l)/n)#a£ + (C7) 

Now, for our purposes, i r = j t and i s = j u since the two states considered here are intended to be the same state 
either translated or acted on by our Hamiltonian (or both) and both of these operations preserve the value of the 
free index. Now, in order to avoid mixing with the singlet states, we also assume i r ^ i s and j t =/= j u . Thus, after 
expanding and contracting indices, the above expression reduces to 



[2(n + l)} (x+x '- 2 ^ 2 [((n + l) d - l)/n + ((n + x)L-d-x-x'+2 - i)/ n ] 
In the large L and large n limits (taking L large first as usual), this reduces to 

2(x+x'-2)/2 L-d-(x+x')/2 



(C8) 



(C9) 



We now consider the case in which t < r < u < s. We label the sites and indices as in Figure ||. We define 
x' = u — t, x — s — r and d = u — r — 1. Routine calculations similar to those above demonstrate that the overlap is 
equal to 



[2(n + i)](*+-'-2)/2-^j : [ 2n (n + l)] d/2 [(((n + i)L- x -x'+d+2 _ i)/ n )^ + sf/, 
Using the same conditions on the indices as above, the overlap reduces to 

[2(n + l)] {x+x '- 2)/2 - d [2n(n + l)] d/2 [{{n + i)i-*-*'+d+2 _ 
In the large L and large n limits (taking L large first as usual), the overlap reduces to 

2(x+x'-d-2)/2 n L+d-(x+x')/2 



(CIO) 



(Cll) 



(C12) 



Lastly, we must consider the case where t < r < s < u. We label the indices as in Figure [lj. We define x' = u — t, 
x = s — r and d = u — s. Routine calculations similar to those show above demonstrate that the overlap is equal to 



jt "t u . 



[2(n + l)} (x '- x)/2 SZSl: [2n(n + l)]^/ 2 [(((n + if-*' - l)/n)«£*£ + 8%5, 
Using the same conditions on the indices as above, this reduces to 



(C13) 
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[2(n + l)] (x '~ x)/2 [2n{n + l)] (x ~ 1)/2 [(((n + 1) W - l)/n)} (C14) 
In the large L and large n limits (taking L large first as usual), the overlap reduces to 

2 (x'-l)/2 Tl i+(x-x')/2-l (C15) 

Thus, we see that the overlap < f\s\t, u > is order n h ~ x if r = t, s = u and higher order in 1/n in every other case. 
Although, we have only dealt with the cases in which s > r and u > t, it can easily be seen that the other cases yield 
the same conclusion. 



APPENDIX D: POSITIVITY OF i/) (x) 

We will now prove that ipo(x) is positive for all odd x. Since we know i/)q(x) is an even function of x, we need only 
prove the result for x > 1. In order for ipg to be normalizable, iI)q{x) must converge to zero. Now, by construction, 
■0o(l) > 0. Assume i/jq(x) > \/x such that 1 < x < X, with x, X odd. We will demonstrate that ip(X + 2) > by 
contradiction. Thus, by induction, ipo(x) > Vx odd. 

Assume x/j(X + 2) < 0. Then, by ( p?S| ) 

MX + 4) = [(X + 2 - nSE - l)/(2V2)]^o(A + 2) - i/j (X) 
MX + 4) - MX + 2) = [{X + 1 - 2V2- nSE)/(2V2)]MX + 2) - ^oPO < 

=> ip{X + A) <ip(X + 2) <0 (Dl) 

Since A + 1 - 2a/2 - n<5£ > 2 - 2\/2 - n5E > 0, t/> (^ + 2) < and M x ) > 0- Now, if M% + 2) < ^ (a;) < for 
x > 3, then 

^ (x + 4) = [(x + 2 - n5£ - l)/(2V2)}4> Q (x + 2) - Vo(x) 
=> Vo(x + 4) - Vo(x + 2) = [{x + 1 - 2V2 - nSE)/(2V2)]Mx + 2) - ^> (x) 
Vo(a; + 4) - Vo(a: + 2) < Vo(£ + 2) - ^ (a;) < 
=*> ^(x + 4) < ^(2; + 2) < 

(D2) 

Since .t + 1 - 2y/2 - nSE > 4 - 2%/2 - nSE > 2y/2. Thus, by induction, M x ) < ip (X + 2) < Vx > X + 2 which 
contradicts the fact that tpo(x) — > as a; — > 00. Therefore, we must have that ipo(X + 2) > which completes our 
initial induction and so ipo{x) > for all x. Here we have used the numerically determined value of nSE ss —3.747. 
Higher energy eigenfunctions will not be positive. 
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